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Data from Au + Au interactions at -^sjvjv = 130 GeV, obtained with the PHENIX detector 
at RHIC, are used to investigate local net charge fluctuations among particles produced near mid- 
rapidity. According to recent suggestions, such fluctuations may carry information from the Quark 
Gluon Plasma. This analysis shows that the fluctuations are dominated by a stochastic distribution 
of particles, but are also sensitive to other effects, like global charge conservation and resonance 
decays. 

PACS numbers; 25.75.Dw 



The PHENIX detector at the Relativistic Heavy-Ion 
Collider (RHIC) is a versatile detector designed to study 
the properties of nuclear matter at extreme temperatures 
and energy densities, obtained in central heavy-ion col- 
hsions at ultra-relativistic energies. A central goal of 
these studies is to collect evidence for the existence of the 
Quark-Gluon Plasma (QGP) characterized by deconfined 
quarks and gluons. 

There are several proposed ways to experimentally ver- 
ify the existence of a QGP . A general problem is that 
many of these signals also can be produced in a hadronic 
scenario, albeit special conditions of highly compressed 
matter have to prevail. Furthermore, it is not straight- 
forward to determine how the various plasma signals are 



distorted when the deconfined matter transforms back to 
hadronic matter. Recent theoretical investigations ||^-||] 
predict a drastic decrease of the event-by-event fluctu- 
ations of the net charge in local phase-space regions as 
a signature of the plasma state. These fluctuations are 
not related to the transition itself, but rather with the 
charge distribution in the primordial plasma state. The 
basic idea is that each of the charge carriers in the plasma 
carries less charge than the charge carriers in ordinary 
hadronic matter. The charge will thus be more evenly 
distributed in a plasma. The main concern of the theo- 
retical discussions is how and why the original distribu- 
tion survives the transition back to ordinary matter [^,0 . 
Predictions, for a rapidity coverage Ay > 1, range up to 
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an 80% reduction in the magnitude of the fluctuations, 
as measured by the variance of the net charge. 

Decays of hadronic resonances influence the net charge 
fluctuations, whether or not deconfinement is reached. 
In the absence of a QGP, deviations from statistical be- 
haviour can be used to determine the abundance of e.g. 
p and uj mesons In a hadron gas resonances are ex- 
pected to decrease the fluctuations by about 25%. Glob- 
ally, fluctuations will be further reduced, since charge is 
a conserved quantity. Deviations from a statistical be- 
haviour in the net charge are also of relevance when the 
phase transition is treated as a semiclassical decay of a 
Polyakov loop condensate . Although multiplicity fluc- 
tuations have been studied extensively in both hadronic 
and nuclear processes [|o|, net charge fluctuations have 
not been addressed experimentally. 

In this Letter we report results from an analysis of 
net charge fluctuations for particles produced in Au+Au 
interactions at ^snn — 130 GeV. The fluctuations are 
studied in the variables R = n+ / n_ , the ratio between 
positive and negative particles, and Q = n+ — n-, the net 
charge The advantages and disadvantages of these 
variables will also be discussed. 

Information from one of the PHENIX central track- 
ing arms (west) is used in this analysis, where events are 
required to have a well-defined vertex close to the cen- 
ter of the apparatus {\Z\ < 17 cm), as defined by the 
two beam-beam counters (BBC). These are Cherenkov- 
counters surrounding the beam, placed on both sides 1.44 
m from the interaction region, covering the pseudorapid- 
ity region 3.0 < \r]\ < 3.9. Together with the informa- 
tion from the two zero-degree calorimeters (ZDC), placed 
further away (18 m), the BBC information is used for 
off-line centrality selection A total of about 5x10^ 
minimum bias events has been analyzed. The PHENIX 
west arm spectrometer has an acceptance of 0.7 units of 
pseudorapidity (-0.35 < ry < 0.35) and 7r/2 radians in 
azimuth ip. Charged-particle trajectories are recorded in 
a multiwire focusing drift chamber (DC) jlj]. The com- 
bination of reconstructed DC tracks |lj] with matching 
hits in the innermost pad-chamber plane (PCI) defines 
the sign of the charge of the particle and also provides 
a high resolution measurement of the transverse momen- 
tum pt of tracks originating from the collision vertex. 
Tracks with a reconstructed pr less than 0.2 GeV/c have 
been excluded from the analysis due to a low reconstruc- 
tion efficiency and large contributions from background 
sources, as revealed by simulations. 

The tracking efficiency and the charge assignment have 
been studied using GEANT |l|] simulations. Of partic- 
ular importance in this context is a realistic description 
of the drift chamber response. The drift distances have 
been calculated based on a geometric model of the indi- 
vidual drift cells. Additional parameters describing the 
response of the chamber, i.e. single-wire efficiency, pulse 
width, single-hit resolution, and space drift-time relation, 



have been extracted from measured data, parameterized, 
and applied empirically in the simulation. 

RQMD simulations are used to study the detection 
efficiency, and the fraction of reconstructed particles that 
preserve their charge, as well as to evaluate the results of 
the analysis. The charge fluctuations in RQMD are con- 
sistent with calculations based on other hadronic mod- 
els like UrQMD and HIJING g. The overaU efficiency 
for detecting a charged particle within the acceptance 
is found to be around 80% for both positive and nega- 
tive particles. Depending on px, between 70% and 85% 
of the reconstructed tracks are in one-to-one correspon- 
dence with a primarily produced particle. The remaining 
tracks come from secondary interactions in the detector 
material and from decays, where the original charge in- 
formation is lost. 

In each event the numbers of positively charged parti- 
cles n-|_, negatively charged particles n_, and their sum 
rich are recorded. In a stochastic scenario, with a fixed 
number of charged particles within the acceptance, where 
each particle is assigned a random charge (-f 1 or —1 with 
the same probability), the variance of the net charge, Q, 
is 



V{Q) = {Q^)-{Qf =nch- 
The normalized variance in Q is 

.«3) = ^ ^ 1. 

rich 



(0.1) 



(0.2) 



For the charge ratio, in the stochastic scenario, V{R) = 
{R^) ~ (R) will approach the value A/uch as Uch increases 
and v{R) = Uch ■ asymptotically approaches 4. A 

small asymmetry between positive and negative parti- 
cles affects v{R) drastically, whereas the effect on v{Q) 
is negligible. If we write the probability , that a given 
particle has positive charge, in the form p+ = Y2(l +e), 
and subsequently P- = V2(l — e), we find 



viQ) 
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(0.3) 



while the asymptotic value of v{R) is 4 -|- 16e -|- 0(£^). 
Detector or reconstruction inefficiencies do not influence 
those results in the stochastic scenario. To calculate v{Q) 
as a function of multiplicity, Eq. 0.3 can be used and v{R) 
can be calculated from 
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and 
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where A = 1 — p""*" — p"!"'' is the new normalization due 
to discarding events with n+ or n_ equal to zero. The 



3 



variance of R, even for a purely stochastic charge distri- 
bution, depends on multiplicity and on the fractions of 
positive and negative particles. 

The data show a small excess of positive particles, 
growing proportionally with rich, in qualitative agree- 
ment with calculations using RQMD and GEANT. A 
part of this excess comes from the intrinsic isospin asym- 
metry and a part from secondary interactions in the de- 
tector and surrounding materials. 

In Fig. la, v{R) and v{Q) are displayed as functions 
of rich- v{Q) is multiplied by a factor of 4 to compensate 
for the asymptotical difference between v(R) and v(Q). 
Both v{Q) and v{R) are well described by the results ob- 
tained from the stochastic scenario, including the positive 
excess, as given by the curves. 

Since v{Q) is independent of rich one expects v(Q) to 
be close to unity also in representations where other cen- 
trality measures are used. On the other hand, since v{R) 
depends on multiplicity, it will have a complicated be- 
haviour as a function of centrality, making it difficult to 
draw any further conclusions. We will thus focus on v(Q) 
for the rest of this analysis. 

In Fig. lb, v{Q) is displayed as a function of central- 
ity based on the ZDC/BBC information. The full event 
sample, corresponding to 92% of the inelastic cross sec- 
tion ||ll[, is divided into 20 centrality classes, where each 
class corresponds to 5% of the events. Class 20 repre- 
sents the most central events. With the increased reso- 
lution on the y-axis in Fig. lb, it is evident that v(Q) 
is consistently below unity, and deviates from stochastic 
behaviour. The value is, however, far above the most 
optimistic QGP predictions v{Q) ~ 0.2 ||^, although one 
should keep in mind that our coverage in rapidity is on 
the limit for these predictions and that we have only par- 
tial coverage in azimuth. 

There may be other explanations for deviations from 
stochastic behaviour than the one offered by the quark- 
gluon plasma. These include global charge conservation 
and neutral resonances decaying into correlated pairs of 
one positive and one negative particle. Both of these 
effects will decrease the fluctuations, and the decrease 
will grow in proportion to the experimental acceptance. 
In a stochastic scenario, taking global charge conserva- 
tion into account, the normalized variance v{Q) becomes 
(1 — p), where p is the fraction of observed charged parti- 
cles among all charged particles in the event. Eventually, 
if all charged particles are detected, v{Q) will become 0. 

Experimentally we can change the fraction p of parti- 
cles within the acceptance by using different regions of 
the PHENIX west arm. In Fig. 2, v{Q), for the 10 % 
most central events, is displayed as a function of Aipd, i.e. 
the chosen azimuthal interval of the spectrometer. For 
comparison, the results from RQMD processed through 
GEANT are shown. The data and the simulation show a 
similar trend. Note that the errors given are correlated, 
since the data in one bin are a subset of the data in the 



next. The solid line corresponds to the (1— p) dependence 
discussed above. The linear relationship between p and 
Aifd is estimated from the phase-space distribution of 
particles in RQMD, including effects from reconstruction 
efficiency and background tracks. For larger angles, both 
data and the RQMD results lie consistently below the 
line, which indicates that effects from resonance decays 
are important. 

Due to the influence of the magnetic field the positive 
and negative particles will have different azimuthal ac- 
ceptance. The Atfd study in Fig. 2 thus selects partly 
non-overlapping regions of phase space for positive and 
negative particles. A remedy for this is to use the recon- 
structed iy9-angle for each particle (pr, i.e. the azimuthal 
direction of the particle at the primary vertex, before it 
is defiected by the magnetic field. Figure 3a shows the 
acceptance in transverse momentum and ipr for positive 
and negative particles. By choosing the azimuthal inter- 
val Aifr symmetrically around the center of the accep- 
tance, a better phase space overlap is achieved for small 
azimuthal intervals. In Fig. 3b, v{Q), for the 10 % most 
central events, is displayed as a function of Aipr. The 
{I — p) dependence, which is no longer linear, is given by 
the solid curve. Again data and the RQMD results show 
a similar trend, but the deviations from the curve are 
larger in this representation, indicating that an overlap 
in phase space is of importance. For large values of Aip, 
the acceptance approaches the limits determined by the 
boundaries of the tracking arm, and the two representa- 
tions are essentially the same. 

The effects of the detector inefficiency and background 
tracks not assigned the correct charge have been investi- 
gated in a Monte Carlo simulation. It is assumed that the 
inefficiency independently removes positive and negative 
particles with the same probability, and that the back- 
ground consists of uncorrelated positive and negative par- 
ticles. The reconstruction efficiency and the amount of 
background have been determined from the RQMD and 
GEANT simulations discussed earlier. Both the ineffi- 
ciency and the background contribution have the effect of 
diluting the signal and pushing the value of v{Q) closer 
to 1. The dilution due to these effects can be treated 
as an experimental systematic error, estimated from the 
simulations, setting a lower limit on v{Q). For the net 
charge fluctuations in the region -0.35 < rj < 0.35, pr > 
0.2 GeV/c, Aip = 7r/2, 

v{Q) = 0.965 ± 0.007(stat.) - 0.019(syst.) (0.6) 

is obtained for the 10 % most central events. A linear ex- 
trapolation of this value to full azimuthal coverage gives 
a value of v(Q) in the range 0.78 - 0.86, in qualitative 
agreement with calculations from a hadronic gas. For 
comparison, it would be desirable to have a Monte Carlo 
model for the QGP which exhibits the predicted reduc- 
tion in the charge fluctuations and which could be used 
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to study the sensitivity of the method for Hmited accep- 
tance. 

To summarize, we have shown that the data behave in 
an almost stochastic manner. There are also clear indica- 
tions that effects from liadronic decays are seen; the data 
are in good agreement with RQMD calculations, which 
includes the effects of global charge conservations as well 
as neutral hadronic resonance decays. Furthermore, the 
data show no centrality dependence, which is in contra- 
diction to the expectations from a Quark-Gluon Plasma 
scenario. We have clearly demonstrated that the fluc- 
tuations of the charge ratio v{R) and of the net charge 
v{Q) are well understood in a stochastic model. We, 
however, advise against the usage of the proposed R vari- 
able 1^ , since it unnecessarily complicates the evaluation 
of the fluctuations, and the intrinsic decrease of v{R), as 
a function of centrality, can be mistaken as a 'plasma fin- 
gerprint'. The measured value of v{Q) = 0.965 is far from 
the value predicted for a plasma. Even extrapolating the 
linear trend seen in the data in Fig. 2 to full azimuthal 
coverage, renders values of the fluctuations, which are far 
above the proposed values. With the caveat of our lim- 
ited acceptance in rapidity, these results clearly indicate 
either the absence of a plasma or that the proposed signal 
does not survive the transition back to hadronic matter. 
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FIG. 1. a) The normalized variances v{Q) and v{R) 
as functions of rich, together with curves for stochastic be- 
haviour, b) The normalized variance v{Q) for different cen- 
trality classes, as described in the text. 
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FIG. 2. v{Q), for central events, as a function of the az- 

imuthal coverage of the detector for data and events simulated 
with the RQMD model. The solid line shows the expected re- 
duction in v{Q) in the stochastic scouario when global charge 
conservation is taken into account. (Angles in degrees.) 
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FIG. 3. a) The acceptance in reconstructed transverse mo- 
mentum pt and azimuthal angle ipr. The solid curves indicate 
the acceptance bands for positive and negative particles, re- 
spectively, b) The effect on v{Q) of varying the acceptance in 
ipr ■ The solid curve shows the expectation from global charge 
conservation. The 10 % most central events are used for data 
and RQMD. (Angles in degrees.) 
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